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ABSTRACT 


A fission -spectrum neutron -dose -rate -attenuation curve was computed by the 
discrete-ordinates method for Li^H and natural LiH for penetration depths up to 
139. 5 cm. Results were compared with a recent Monte Carlo calculation for natural 

£ 

LiH. A single value for an energy-independent-neutron-removal cross section for Li 
and natural Li in LiH was deduced for use in an Albert -Welton dose -attenuation kernel 
but is accurate for only a given range of penetration. The validity of fast -neutron- 
removal theory was demonstrated for the cases of W followed by LiH. From this pene 
tration data, a neutron -removal cross section for W was calculated. 



FISSION NEUTRON ATTENUATION IN LITHIUM-6, NATURAL 
LITHIUM HYDRIDE, AND TUNGSTEN 
by Gerald P. Lahti 
Lewis Research Center 

SUMMARY 

The calculated attenuation of dose due to a point source of uranium -2 35 fission- 
spectrum neutrons through lithium -6 hydride (Li H) is presented. The computations 
were made with the use of a one -dimensional discrete -ordinates transport code. The 
single -scatter rad dose was calculated to a depth of 139. 5 centimeters. A similar cal- 
culation was performed for natural lithium hydride (LiH) and compared with recently 
published Monte Carlo results for the same configuration. The results are within 
25 percent of the smoothed Monte Carlo results at deep penetration depths. Li H and 
LiH offer similar attenuation properties; because of its 12-percent lower density, Li H 
offers a reduction in weight over LiH. 

An attempt was made to deduce a single energy -independent neutron-removal cross 

section for Li and Li in LiH for use in an Albert -Welton dose -attenuation kernel utilized 

by point-kernel shielding codes. Although no single value gave close agreement to the 

calculated results over the entire range, values of 0. 100 and 0. 115 square centimeter 

R 

per gram for natural Li and Li , respectively, in LiH gave the best agreement. A 
better fit of the dose attenuation is afforded by a dose -attenuation kernel for LiH of the 
form a Q exp(ajr + a 2 r + a 3 r ), where the a's are constants and r is the distance 
from the source. 

The validity of the fast-neutron-removal theory for cases of tungsten (W) followed 
by Li H or LiH is demonstrated for cases of a poipt -fission source surrounded by up to 
25. 5 centimeters of W followed by Li H or LiH. A fast -neutron-removal cross section 
of 0. 0105±0. 001 square centimeter per gram for W provides a good fit to the calculated 
data. 


INTRODUCTION 


Lithium -6 hydride (Li H) is a neutron shield material of interest because the use of 


the lithium-6 isotope with its high (n,a) cross section results in the depression of 
secondary gammas (arising from (n,y) events) in shields involving both neutron and 
gamma shield arrangements. Also, Li®H has a lower density than that of natural 
lithium hydride (LiH), so that it offers the additional possibility of reduction in shield 
weight, provided that it is as good a fast -neutron (dose) attenuator as LiH. 

A series of discrete-ordinates calculations was made to compare the fast -neutron- 
dose attenuation of LiH and Li H and to generate removal cross sections for use in 
point -kernel shielding codes (e. g. , QAD, refs. 1 and 2). The calculations were made 
with the use of the one -dimensional discrete -ordinates transport code ANISN (ref. 3) 
for a small (point) source emitting a fission spectrum of neutrons in an essentially in- 
finite medium of LiH and Li H. A comparison was made between the discrete-ordinates 
calculation for LiH and a recent Monte Carlo calculation (ref. 4) made for the same 
configuration. Similar calculations were made for the cases of tungsten (W) followed by 
Li H or LiH to demonstrate the validity of fast -neutron -removal theory and to obtain a 
fast -neutron -removal cross section for W. 


FISSION NEUTRON ATTENUATION IN LITHIUM -6 HYDRIDE AND LITHIUM HYDRIDE 

Details of Input to Calculations 

The discrete -ordinates code ANISN was used to make transport calculations for the 

235 

cases of a uranium-235 (U ) fission-spectrum source located at the center of a 139. 5- 
centimeter-radius sphere of LiH or Li 0 H. The computations were made with the use of 
multigroup cross sections generated by GAM n (ref. 5) for the 25 fast energy groups in- 
dicated in table I. 

The energy group structure is designed to place emphasis on that energy region con- 
tributing most to the neutron dose. A more detailed discussion of the energy group 
structure may be found in the section, Importance of Fine Group Structure Above 6 MeV. 
A single thermal group cross section was added to complete the considered energy span 
from 14. 9 to 0. 0 MeV. A 16 -point full-range Gauss -Legendre quadrature scheme was 
used to describe the angular fluxes. The Pg approximation was used to describe the 
elastic -scatter-source term. The mesh spacing used is indicated in table H. 

The v 0 fission-spectrum source (ref. 5) was assumed to be spread over the first 
five mesh intervals to facilitate the computation. The flux-to-dose conversion factors 
used were taken from figure 1 of reference 6 and are listed in table I. The resulting 
dose is a first-collision-rad dose. The density of natural LiH used was 0. 75 gram per 
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TABLE I. - ENERGY GROUP STRUCTURE AND FLUX-TO- 


DOSE CONVERSION FACTORS 


Group 

Lower energy limit 

Flux-to-dose conversion 


of group, 

factor, 4jrr 2 xDose, 


E, 

(rad/hr)/ (neutron/ (cm 2 ) (sec)) 


eV 


1 

13. 5X10 6 

2. 05x10”® 

2 

12.2 

2.05 

3 

11. 1 

2.05 

4 

10.0 

2.05 

5 

9.0 

2.05 

6 

8. 2 

2.05 

7 

7.4 

1.84 

8 

6.7 

1. 84 

9 

6. 1 

1.84 

10 

5. 5 

1.76 

11 

5.0 

1.69 

12 

4. 5 

1.66 

13 

4. 1 

1.62 

14 

3.7 

1.62 

15 

3.0 

1.55 

16 

2.5 

1. 37 

17 

2.2 

1. 15 

18 

1.8 

1.12 

19 

1. 35 

1.08 

20 

1. 10 

.935 

21 

.91 

. 83 

22 

. 55 

.72 

23 

. 11 

. 54 

24 

3.0X10 3 

.072 

25 

. 4X10° 

.0036 

26 

0 

.004 


TABLE H. - MESH SPACING FOR DISCRETE- 
ORDINATES PROBLEM 


Range of penetration 

Mesh interval 

depth, 

size, 

r, cm 

cm 

0 to 6.0 

0.4 

6.0 to 10. 5 

.75 

10. 5 to 139. 5 

1.00 


3 



TABLE m. - DENSITIES OF LITHIUM -6, 
LITHIUM-7, AND HYDROGEN USED IN 


LITHIUM HYDRIDE CALCULATIONS 


Constituent 

Density, p 

isotopes 
in - 

g/cm 3 

atoms/cm 3 

Natural lithium 
hydride: 
Hydrogen 

0.0951 

0.0 56837X10 24 

Lithium -6 

.0421 

.004217 

Lithium -7 

.6130 

.052619 

Total 

0.750 


Lithium -6 hydride: 


0. 056837X10 24 

Hydrogen 

0.0951 

Lithium -6 

.5677 

.056837 

Total 

0.6628 



cubic centimeter. The density of each of the constituent isotopes making up LiH and 

Li°H is listed in table m. The density of Li H was calculated to be 0. 6628 gram per 

6 

cubic centimeter. The hydrogen atom density in Li H was assumed to be unchanged from 

7 fi 

that in the natural LiH, and the Li atom was assumed to be replaced by Li. 


Results of Discrete-Ordinate Calculations for Lithium Hydride 

Comparison of attenuation propert ies of lithium^ hydride and lithium -6 hydride. - 

Sernilogarithmic plots of the attenuation of neutron dose in LiH and Li b H are shown in 

2 2 

figure 1. The expression 4nr x Dose ((cm ) (rad/hr)/ (fission neutron/sec)) is plotted 
against r, the distance from the center of the source. (All symbols are defined in 
appendix A. ) 

The dose in Li H is about 10 percent below that of LiH for r of about 30 centi- 
meters, is about the same at 65 centimeters, and is about 60 percent above that of LiH 
at a depth of 120 centimeters. Small differences in neutron attenuation by Li H and LiH 
have been reported for experiments performed with 4- and 12-inch (10. 16- and 30. 48-cm) 
(ref. 7) and 36-inch (91.44-cm) (ref. 8) slabs of LiH and Li^H. Both curves in figure 1 
are concave upward after a penetration of about 30 centimeters, which indicates a con- 
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Figure 1. - Neutron dose attenuation in lithium hydride as calculated by discrete- 
ordinates method. 


tinually hardening spectrum. Table IV lists the local attenuation coefficient p. where 

_o 

the dose attenuation of figure 1 is assumed to be of the form Ar exp(-/ir), or in terms 
of the density p and the mass thickness pr, Ar - ^ exp [-(p/p) (Pr)], where A is a con- 
stant. The local attenuation coefficient decreases with increasing depth in LiH. At a 
depth of 130 centimeters, each curve of figure 1 appears to attain a constant slope de- 
fined by a single exponential. The value of this slope, as indicated in table IV, is 
0. 112 per centimeter (0. 168 cm^/g) for Li^H and 0. 1115 per centimeter (0. 149 cm^/g) 
for LiH. 

Shielding effectiveness of lithium-6 hydride and lithium hydride . - With respect to 
fast-neutron-dose attenuation, Li®H is favored over LiH because of the slight weight 
advantage it offers. For example, both afford the same attenuation at a 60-centimeter 
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TABLE IV. - VALUES OF LOCAL DOSE -ATTENUATION COEFFICIENTS 




0 20 40 60 80 100 

Penetration depth in lithium hydride, r, cm 


Figure 2. - Neutron dose attenuation in lithium hydride 
as calculated by Monte Carlo method. Lithium hydride 
density, 0.75 gram per cubic centimeter. 
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depth, but Li H is 88 percent as dense as LiH (a 12 -percent weight reduction). However, 
the weight advantage may decrease at deep penetration depths. 

Com parison of discrete-ordinate results wi t h Monte Carlo results. - Shown in fig- 
ure 2 is a semilogarithmic plot of the expression 4irr% x Dose ((cm^) (rad/hr)/ 

(fission neutron/sec)) against r for LiH as calculated by ANISN. Also included are the 
Monte Carlo results of Kam and Clark (ref . 4) shown with 3 -standard -deviation error 
bars and a smoothed curve of the Monte Carlo results. The discrete -ordinates results 
are virtually the same as the smoothed Monte Carlo results to a depth of 40 centimeters. 

For deeper penetrations, the present calculations are lower than the smoothed Monte 
Carlo results by, at most, 25 percent. 

A calculation was made to ascertain the effect of the approximation to the elastic - 
scatter-source term on the discrete-ordinate -calculation results. The approxima- 
tion result (for both Li 6 H and LiH) is 10 percent below the P 3 result at 40 centimeters 
(2j decades of attenuation) and 50 percent lower at 120 centimerers (about 7^ decades of 
attenuation). To obtain accurate results at deep penetration depths requires the P 3 
scattering approximation; at shallow penetrations, P^ scattering is adequate. 

Another possible source of discrepancy in these results is due to the choice of the 
multigroup structure, which is discussed in more detail in the section Importance of 
Fine Group Structure Above 6 MeV. 


Curve Fits to Discrete Ordinates Results 

Evaluati on of lithium removal cross section . - A dose kernel that is often used in 
point -kernel shielding codes to estimate the fast-neutron dose is the Albert -Welton 
kernel, which is of the form 

a. 

4yrr 2 D = e‘ c h a 2 ex P(-“ 3 h ) 


where 


D 

a v a 2 , 

a v a 4 


dose rate from a unit source of fission -spectrum neutrons at a distance r 
coefficients used to account for attenuation in hydrogen 



7 


removal term accounting for attenuation due to all nonhydrogenous material 

O 

equivalent thickness of hydrogen with density of 0. Ill g/cm encountered 
between source and receiver point 

xi. 

fast-removal cross section for M element 

xi. 

mass thickness of M in element 

One set of coefficients a that describe neutron penetration in hydrogen are 
Q!^ = 8.35x10”^ ((cm 2 ) (rad/hr)/ (neutron/sec)) 

a 2 = 0. 29 (dimensionless) 

Og = 0. 83 (dimensionless) 

a ^ = 0. 58 (dimensionless) 

A value of 2 R /p for lithium is deduced for use in the Albert -Welton kernel by 
assuming the framework of this kernel to be valid for LiH with 2p/p as the only ad- 
justable parameter. (Refs. 7 and 8 have also estimated the value of 2 R /p for Li in 
LiH with the use of Albert -Welton kernel as a framework for a curve fit. ) Figure 3 is 
plot of Albert-Welton kernel evaluations for several values of S R /p. Also shown for 
reference is the appropriate Pg discrete -ordinates result. Figure 3(a) illustrates 
Albert-Welton kernel evaluations for values of 2 R /p = 0.0876 (refs. 9 and 10), 0. 100 
(ref. 7), and 0. 107 (ref. 8) square centimeter per gram for LiH*. The value of 0. 100 
square centimeter per gram provides the best fit to the discrete -ordinates results for 

f* 

natural LiH. For the matching of the Li H Pg calculation (fig. 3(b)), a value of 
S_/p = 0. 115 square centimeter per gram fits well from penetration depths of 40 
to 90 centimeters. Accurate fits are obtained in either LiH or Li H over a limited 
range of penetrations. 



*The 2 p/p of refs. 9 and 10 is a slab-removal cross section determined in early 
Lid -Tank experiments. Also, the value of 2 R /p of 0. 100 cm 2 /g was deduced from 

ref. 7 with the assumption that a density of 0. 75 g/cm 2 for LiH (Li density is 0. 655 

3 1 

g/cm ) is to be used with the reported value of 2 D = 0. 0655 cm . Similarly, ref. 8 

1 ** 99 3 

indicates a S R s 0. 07 cm for Li in LiH with a density of 5. 64x10 atoms/cm . 
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Fast-removal cross section, 

Zr'p. 


0 

0. 0876 ^ . 

□ 

100 Albert-Welton 

A 

* 107 > kernel 

O 

* 115 evaluations 

O 

! 120 , 


■ P3 Transport calculation 



40 6 0 80 100 120 0 20 40 60 80 100 

Penetration depth in lithium hydride, r, cm Penetration depth in lithium -6 hydride, r, cm 

(a) Natural lithium hydride. (b) Lithium-6 hydride. 

Figure 3 . - Curve fits for neutron attenuation in lithium hydride. 


TABLE V. - COEFFICIENTS FOR LITHIUM HYDRIDE 
DOSE KERNEL (EQ. (1)) 


Coefficient Lithium hydride Lithium -6 hydride 


a 0 4.1595X10' 5 

3.6550X10" 5 

aj -1.8923X10' 1 

-1. 9597x10 1 

a 2 3.2378X10' 4 

5.1652X10' 4 

a Q 3.802£<10' 8 

-7. 8209X10' 7 



Alternate fit to lithium hydride attenuation_curve. - As shown in figure 3, the 
Albert -Welton kernel does not fit either the magnitude or the slope of the calculated dose 
curve over the entire interval. A better fit to the curve of LiH attenuation for both 
magnitude and slope is offered by a function of the form 

2 2 3 

47rr x Dose = a.Q exp(a^r + a 2 r + a 3 r ) (1) 

where r is the depth (cm) of LiH. Coefficients for LiH and Li H were calculated by 
use of a least -squares fit to the data curves of figure 1 and are listed in table V. 

Equation (1), with the appropriate coefficients, obtains values of 4nr^ x Dose which 
are within ±5 percent of the values calculated with the transport code over the range of 
r from 20 to 139. 5 centimeters (fig. 1). An attenuation kernel of the form of equation 
(1) could be employed in a point -kernel code for layered configurations of high Z -layers 
with LiH. The addition of a removal term e” c as described previously for the Albert - 
Welton kernel, is required to account for fast -neutron attenuation in high Z -layers. 


FISSION NEUTRON ATTENUATION IN TUNGSTEN AND LITHIUM HYDRIDE 

The neutron-removal theory used in neutron shielding calculations may be stated as 
the following equation for the case of a point source of fission-spectrum neutrons sur- 
rounded by a high Z- (atomic number) material followed by a hydrogenous material: 

47rrj x Dose (t z ,t H ) = 47rr 2 X Dose (0,t H ) exp ^ (pt) z (2) 


where 


r l 


” t H + t Z 


Dose (t z >tjj) dose observed at distance r^ from point source with t z cm of high 
Z -material followed by tjj cm of hydrogenous material surrounding 
the source 


H 


thickness of high Z -material, cm 
thickness of hydrogenous material, cm 
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TABLE VI. - ALTERNATE ENERGY GROUP STRUCTURE 


AND FLUX-TO-DOSE CONVERSION FACTORS 


Upper energy limit 
of group, 

E, 

eV 

Flux -to -dose conversion 
o 

factor, 4ttt x Dose, 

2 

(rad/hr)/ (neutron/ (cm )(sec)) 

14. 92X10 6 

2.05X10' 5 

12.21 

2.05 

10.00 

2.05 

8.19 

1.84 

6.07 

1.76 

4.97 

1.66 

4.07 

1.51 

3.01 

1. 37 

2.47 

1.15 

2.23 

1.12 

1.83 

1.08 

1.35 

.94 

1.11 

.83 

.91 

.72 

.55 

.65 

.41 

. 54 

.H 

• 17 

1. 50X10 4 

3.2X10' 7 

3. 35X10 3 

7. 2xl0~ 8 

5. 83X10 2 

1. 3xl0‘ 8 

1.01X10 2 

3.6X10' 9 

2.90X10 1 

3. 6x10 ‘ 9 

1.07X10 1 

5. 4x10 ' 9 

3.06X10° 

9. 5X10 ‘ 9 

1.13 

1.5X10' 8 

.41 

4. OxlO -8 



Dose (0,t H ) dose observed at distance rg from point source with t R cm of hydro- 
genous shield surroimding only the point source 

2 

(St ,/ p) high Z -material-effective removal cross section, cm /g 

±t z 

O 

p density of high Z -material, g/cm 

(pt) z = p x t z 

This formulation simply separates the dose -attenuation kernel into two parts: one 
is due to the high Z -material, and the other is due to the hydrogenous material, which 
functions as if the high Z -component were not present. (Details of this formulation may 
be found in refs. 11 or 12.) This theory has been investigated by discrete -ordinates 
calculations for cases of several thicknesses (6, 10. 5, 20. 5, and 25. 5 cm) of W at a 

c 

density of 19. 3 grams per cubic centimeter followed by Li H or LiH. In these calcula- 
tions, the mesh spacing indicated in table n and the energy group structure of table VI 
were used. The energy group structure is revised from that of table I so that neutrons 
slowing down from the inelastic collisions with W might be followed more accurately. 

At deep penetrations, the difference in neutron dose resulting from the use of this choice 
of energy group structure instead of that of table I is less than 10 percent at a depth of 
90 centimeters for the case of figure 1. 

Results of Discrete-Ordinates Calculations for Tungsten and Lithium Hydride 

Effective attenuation of fission neutrons by tungsten. - Figure 4 is a plot of 
47rr 2 X Dose ((cm 2 ) (rad/hr)/ (fission neutron/sec)) against the depth of penetration in the 
LiH following W, as obtained by the discrete-ordinates calculation. Figure 4(a) is for 
W thicknesses of 0, 6, 10.5, 15.5, 20.5, and 25. 5 centimeters followed by Li H. Fig- 
ure 4(b) is for W thicknesses of 0, 6, 10. 5, and 15. 5 centimeters -followed by natural 
LiH. If the removal theory formulation is valid, the resulting curves should be sepa- 
rated from the case of LiH without the presence of W by the constant factor 
exp [(S R /p) (pt) w J where (S R /p) is the removal cross section for W and (pt)^ is 

the mass thickness of W. Inspection of figure 4 shows that an approximately constant 
factor separates the two curves after a LiH penetration of about 20 centimeters, a min- 
imum penetration required to restore the spectrum to that characteristic of LiH. 

Table VH lists ratios of dose without W to dose with W (the apparent attenuation for a 
given thickness of W) for several depths in LiH as calculated by the discrete -ordinates 
code. Also included in the table is the effective removal cross section for W that is 
calculated from equation (2) . 
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Neutron dose attenuation, 47 rr 2 xDose, (cm2)(rad/hr)/(neutron/sec) 




(a) Lithium-6 hydride. <W Lithium hydride. 

Figure 4. - Discrete-ordinates calculation of fission neutrons in tungsten followed by lithium hydride. 



TABLE VH. - APPARENT ATTENUATION AND EFFECTIVE REMOVAL CROSS SECTION 


FOR TUNGSTEN FOLLOWED BY LITHIUM HYDRIDE 


Penetration 

[Itatio of dose without tungsten 

1 Calculated effective removal cross section, 

depth, 

to dose with tungsten 

a 


S R /p, cmVg 


r, 

o 

4zrr X Dose (0 

’ 4 LiH 







cm 

- 




— 







4jrr X Dose (t^y, ^LiH^ 







Thickness of tungsten, cm 

Thickness of tungsten, t^, cm 

1 


6 

10. 5 

15.5 

20. 5 

25.5 

6 

10. 5 

15.5 

20.5 

1 25. 5 





Lithium - 

6 hydride 





30 

3.41 

9.28 

29.0 

91.8 

292 

0.01060 

0.01100 

0.01126 

0.01140 

0.01153 

40 

3.41 

9.24 

28. 5 

89.0 

280 

.01060 

. 01096 

.01120 

.01134 

.01145 

50 

3.35 

8.91 

27.0 

82.7 

254 

.01045 

.01080 

.01102 

.01115 

.01125 

60 

3. 28 

8. 59 

25.4 

76. 2 

231 

.01025 

.01058 

. 01082 

.01095 

.01106 

70 

3.30 

8.19 

23.8 

66.6 

197 

. 01030 

. 01036 

. 01060 

.01075 

. 01086 

80 

3.16 

7.93 

22.85 

64.2 

189 

.00993 

.01022 

. 01046 

.01060 

.01073 

90 

3.09 

7.76 

22.15 

62.8 

186 

.00973 

.01012 

.01035 

.01050 

.01065 

100 

3.07 

7.64 

21.7 

62.3 

— 

. 00968 

.01003 

. 01029 

.01045 

.01062 

110 

3.05 

7. 58 

21. 5 

— 


. 00960 

. 00999 

.01025 

.01043 







Lithium hydride 





30 

3.43 

9. 31 

29. 1 


— 

0.01062 

0.01101 

0.01126 



40 

3.42 

9. 30 

28.9 

— 

— 

.01058 

.01100 

.01124 



60 

3.24 

8.55 

25.8 

— 

— 

.01015 

. 01060 

. 01086 



80 

3.06 

7.85 

22.8 

— 

— 

. 00966 

.01017 

. 01045 



100 

2.94 

7.42 

21.2 

— 

— - 

. 00931 

. 00990 

. 01020 




a From fig. 4. 


Effective removal cross section for tungsten. - From the foregoing remarks and the 
data of table VII, it may be concluded that simple removal theory will give reasonable 
results for deep penetration problems involving up to 25. 5 centimeters of W and 
110 centimeters of LiH for single W -layer configurations for depths in LiH greater 
than 30 centimeters. The calculated effective removal cross sections for W listed in 
table VII can be compared with the two experimentally observed values of 0.01145± 

0. 00065 and 0. 0102 5±0. 00082 square centimeter per gram and with the weighted average 
values of 0. 0100 square centimeters per gram (given in ref. 10). All the values of cal- 
culated for W fall within the experimental errors reported. An increase occurs 

in the effective removal cross section of W with W thickness and a decrease occurs 
with LiH thickness. The value of Sj^/p = 0. 0105±0. 001 square centimeter per gram 
covers the entire range calculated. 
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Neutron spectra i n tungsten and lithium hydride . - Another interesting result is 

illustrated in figure 5, a plot of the relative integral neutron spectrum (i. e. , the fraction 

o 

of neutron flux above some energy E) for several values of E against the distance 

from the source r. The dashed curves are for the case of a point -fission source in 

Li H, and the solid curves are for the case of a point- fission source surrounded by 

15.5 centimeters of W followed by 124 centimeters of Li H. 

For the case of 15. 5 centimeters of W followed by Li H, the neutron spectrum is 

shifted by nonelastic scattering events in W. The energy spectrum of neutrons leaking 

from the W (at r = 15. 5 cm) is about 40 percent in the energy interval of 110 to 410 KeV 

and another 35 percent in the interval of 15 to 110 KeV, which is much different from the 

fission spectrum. After about 25 centimeters of penetration in Li H, the relative 

£ 

spectra (for the case of W followed by Li H) are virtually the same as for the case of 
/»/»/» 

Li°H alone. (If the solid curve (W-Li°H) were plotted against depth in Li H, it would be 

shifted 15.5 centimeters to the left, which would result in almost the same values that 

are plotted for the dashed (Li H only) curve after the spectrum has reached equilibrium.) 

Regardless of the low-energy-leakage spectrum, the high energy flux, greater than 6 

or 8 MeV, leaking from the high Z -shield (or source) gives rise to the characteristic 

spectrum shape and magnitude at deep penetrations in hydrogenous materials. A tabula- 

£! 

tion of multigroup fluxes at various penetrations of LiH and Li H appears in appendix B. 


Importance of Fine Group Structure Above 6 MeV 

Figure 5 illustrates how the fraction of neutrons above 6 MeV (and the fraction of 
dose attributable to them) increases with penetration depth. At one point in this series 
of computations, an alternate, coarser 14-group structure was used to ascertain its effect 
on calculated dose. This structure had only one energy group in the region extending 
from 6.1 to 14.9 MeV, as compared with the 26 -group set of table VI with four energy 
groups in the region. The dose in the 14 -group result was lower than that for the 
26-group set at the following depths in LiH: 10 percent at 40 centimeters, 45 percent at 
70 centimeters, a factor of 2 at 90 centimeters, and a factor of 8. 5 at 120 centimeters. 
The increasing error in the 14 -group calculation is directly attributable to the lack of 
definition above 6 MeV; more precisely, the calculation lacks a correct description of 
neutron transport characteristics above 6 MeV, that is, neutrons whose energy spectrum 
changes greatly with the penetration depth and is no longer that calculated by the infinite - 
medium cross-sectipn code (GAM-II) for a single group above 6 MeV. At a 120- 
centimeter penetration, each of the four groups above 6 MeV in the 26 -group calculation 


^This fraction is equal to J g°° cp( r, E) dE / -fi cp( r, E) dE. 
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Relative flux greater than energy E 



Figure 5. - Relative flux in tungsten and lithium-6 hydride. 




contributes about 10 percent to the total dose. The average energy of those neutrons 
above 6 MeV (see fig. 5) shifts from 7.48 MeV, as calculated by GAM-H for the infinite 

medium, to 7. 60 MeV at 20 centimeters, to about 8. 5 MeV at 60 centimeters, and to 
9.2 MeV at 100 centimeters. 


CONCLUDING REMARKS 

A fission -spectrum neutron-dose -rate -attenuation curve was computed by the 

fi 

discrete -ordinates method for lithium -6 hydride (Li H) and natural lithium hydride (LiH) 

at penetrations up to 139. 5 centimeters for the case of a point source of fission- 

spectrum neutrons in an essentially infinite medium. Results were compared with a 

Monte Carlo calculation for natural LiH and were in general agreement. On the basic of 

penetration depth, Li H and LiH offer comparable attenuation properties. 

If the Albert -Welton kernel for fast -neutron-dose attenuation is used for LiH or 
fi R ft 

Li H at deep penetrations (10° to 10 attenuation), a value of 0. 100 or 0. 115 square 
centimeter per gram should be used for the fast-removal cross section of natural Li or 

ft 

Li , respectively, in LiH. These values result in an approximate fit to the discrete - 

ordinates results for a given range of attenuation. A better fit is offered by a single 

2 3 

kernel for LiH of the form aQ exp(a^r + + a^r ). 

The removal theory was demonstrated to describe the attenuation of fission -neutron 
dose for cases of a fission -spectrum source surrounded by tungsten followed by Li H 
or LiH after a penetration of about 30 centimeters in LiH. A good fit was provided by 
the description of the tungsten dose attenuation given by the exponential 

exp jj-(Sj^/p)(pt) J with = 0.0105±0.001 square centimeter per gram. This fit is in 

agreement with experimental results. 

Accurate discrete -ordinates results were obtained with scattering and a coarse 
group structure to a depth of 40 centimeters in LiH; for deeper penetration, Pg scatter- 
ing and a fine high-energy (greater than 6 MeV) group structure is required. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 25, 1968, 

124-09-01-03-22. 
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APPENDIX A 


SYMBOLS 


A 

a Q , a^, &2 j 3. ^ 

c 

E 

h 

r 

t 

Z 

a v a 2 ,a 3 ,a 4 

M 

M/P 

P 

Pt 

S R 

V p 

<P 


constant in dose -kernel expression 
constants in alternate dose kernel 

e~ c accounts for neutron attenuation in high Z -material 
energy, MeV 

equivalent hydrogen thickness, cm 

radius (penetration depth or distance from center of source), 
thickness, cm 
atomic number 

constants in Albert -Welton dose kernel 
local dose -attenuation coefficient, cm"* 

9 

local mass dose -attenuation coefficient, cm /g 

q 

density, g/cm 


-1 


mass thickness, g/cin 
fast -neutron -removal cross section, cm 

2 

mass fast -neutron-removal cross section, cm /g 

o 

neutron flux, neutrons/ (cm ^sec) 


Subscripts (used with S R /p, pt , and t): 
th 


M 

W 

Z 


M element other than hydrogen 
tungsten 

high-Z material 


cm 
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APPENDIX B 


NEUTRON FLUX IN LITHIUM HYDRIDE 

Of possible interest are the numerical values of the neutron flux as a function of the 
distance from the source for the case of a point source of fission neutrons (1 fission 
neutron/sec) in LiH, as calculated by the discrete -ordinates method. This appendix con- 
tains values of the multigroup fluxes for values of r of 15, 30, 60, 90, and 120 centi- 
meters and for the group structure of table VI. The values given represent the multi- 
group fluxes 



<p(E, r) dE s <p.(r) 


The values of co, (r) are listed in table VTII(a) for the case of a point source in LiH and in 

1 g 

table vm(b) for the case of a point source in Li°H. 
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TABLE VIII. - MULTIGROUP FLUX AS CALCULATED BY DISCRETE -ORDINATES METHOD 


(a) In lithium hydride. (b) In lithium -6 hydride. 


Group 

Neutron flux, ^(r) 

[neutrons/ (cm 2 ) (sec)] / (source neutron/ sec) 

Group 

Neutron flux, <p.(r) 

[neutrons/(cm 2 ) (sec)] / (source neutron/sec) 


Penetration depth, 

r, cm 



Penetration depth, i 

, cm 


15 

30 

60 

90 

120 

15 

30 

60 

90 

120 

1 

1.47X10" 8 

9. 16X10 -10 

1.36X10' 11 

3.34X10' 13 

9.69X10" 15 

1 

l.saao " 8 

9.77X10' 10 

1. 56X10' 11 

4.11X10' 13 

1. 29x10 ' 14 

2 

6.96X10' 8 

3. 58X10" 9 

3. 70X10' 11 

6.64X10* 13 

1.5fc<10" 14 

2 

7.6^<10' 8 

4. 28X10' 9 

5. 2^10‘ U 

1.08X10' 12 

2. 78 

3 

2. 21xl0‘ 7 

9.17X10' 9 

6.36X10' 11 

8. 36X10' 13 

1.57 

3 

2.40X10" 7 

1.09X10' 8 

9.16X10' 11 

1.43 

3.05 

4 

9. 70X10 ‘ 7 

3. 12<10~ 8 

1. 33X10’ 10 

1.25X10" 12 

2.01 

4 

9. 77X10' 7 

3. 28 

1. 63X10 * 10 

1.86 

3.50 

5 

1.36X10" 6 

3. 74X10' 8 

1.20 

9.37X10" 13 

1.38 

5 

1.27X10' 6 

3.46 

1.25 

1.21 

2. 14 

6 

2.14 

5. 2SX10 -8 

1.38 

9.46X10" 13 

1.32 

6 

1.91 

4. 58X10’ 8 

1.32 

1.15 

1.98 

7 

4.86 

1.08X10 -7 

2.31 

1.42X10" 12 

1.89 

7 

4. 35 

9.46 

2.18 

1.67 

2. 77 

8 

4.42 

9. 39X10" 8 

1.77 

9.88X10' 13 

1.26 

8 

4.18 

8.91 

1. 77 

1.20 

1.90 

9 

2.49 

5. 19X10' 8 

9.22X10' 11 

4.94X10" 13 

6. 18X10” 48 

9 

2.46 

5.25 

9.79X10' 11 

6.20X10* 13 

9.5»<10‘ 15 

10 

5.48 

1. 13X10' 7 

1.90X10' 10 

9.79X10 -13 

1.20X10 -14 

10 

5.40 

1.15X10' 7 

2.05X10' 10 

1. 2&<10" 12 

1. 83x10 ' 14 

11 

9.03 

1.86 

2.93 

1.41X10' 12 

1.68X10" 14 

11 

8.63 

1.85 

3.07 

1. 71X10* 12 

2. 50 

12 

6.19 

1.28 

1.94 

8.91X10' 13 

1.03X10' 14 

12 

5.62 

1.20 

1.92 

1.0^<10* 12 

1.46 

13 

6.07 

1.27 

1.87 

8. 33X10" 13 

9.46X10 -15 

13 

5.23 

1.11 

1.73 

9.07X10' 13 

1.28 

14 

1.39X10’ 5 

2.92 

4.17 

1. 80X10 ^ 2 

2. OlxlO -14 

14 

1.08X10' 5 

2.25 

3.46 

1.79X10 -12 

2.51 

15 

6. 86X10 “ 6 

1.44 

2.03 

8. 66X10” 13 

9.59X10" 15 

15 

4.96X10' 6 

1.02 

1.56 

8.03X10' 13 

1.13 

16 

2.02X10’ 5 

4.17 

5.88 

2.51X10 -12 

2.78X10” 14 

16 

1.23X10* 5 

2.50 

3.86 

2. OlxlO -42 

2.83 

17 

2.07 

4.24 

5.87 

2.47 

2.70 

17 

1.23X10' 5 

2.46 

3.74 

1.91X10' 12 

2.67 

18 

1.32 

2.68 

3.68 

1.52 

1.65 

18 

6.69X10' 6 

1.33 

2.00 

1. 01X10 -12 

1.40 

19 

1. 54 

3.10 

4.22 

1.73 

1.86 

19 

5.68X10' 6 

1.12 

1.67 

8.42X10' 13 

1.16 

20 

1.50 

3.01 

4.05 

1.64 

1.75 

20 

3.03X10' 6 

5.95X10' 8 

8.88X10 -11 

4. 45X10 ' 13 

6.15X10' 15 

21 

8.87X10" 6 

1.79 

2.64 

1.06 

1.13 

21 

9.21X10' 7 

1.81X10' 8 

2.69X10" 11 

1.35X10' 13 

1. 86X10" 15 

22 

6.92 

1. 38 

1.84 

7. 36X10" 13 

7. 78X10" 15 

22 

2.92X10' 7 

5.73X10" 9 

8. 53X10' 12 

4. 29X10' 14 

5. 90X10' 16 

23 

6. 79 

1.35 

1.79 

7.16 

7.55 

. 23 

1.15X10' 7 

2. 26X10 -9 

3. 36X10" 12 

1.69X10' 14 

3. 32xl0" 16 

24 

3.71 

7. 38xl0“ 8 

9.77X10“ U 

3.89 

4.10 

24 

2.45X10" 8 

4.8^<10' 10 

7.17X10' 13 

3. 59X10' 15 

4.95X10" 17 

25 

2. 24 

4.46X10' 8 

5.90X10" 11 

2. 35 

2.47 

25 

6. 75X10 ‘ 9 

1.33X10' 10 

1.98X10' 13 

9.92X10' 16 

1.37X10 -17 

26 

2.43 

4.83X10" 8 

f 6.40X10" 11 

2.54 

2.68 

26 

2. 39X10 " 9 

4. 70X10' 11 

6.97X10' 14 

3.52X10' 16 

4.81X10" 18 
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